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Abstract 

We constrain the parameters of the variable Chaplygin gas model, using the lo- 
cation of peaks of the CMBR spectrum and SNe la "gold " data set. Equation of 
state of the model is P = —A(a)/p, where A(a) = Aoa~ n is a positive function of the 
cosmological scale factor a, Aq and n being constants. The variable Chaplygin gas 
interpolates from dust-dominated era to quintessence dominated era. The model is 
found to be compatible with current type la Supernovae data and location of first peak 
if the values of Sl m and n lie in the interval [0.017, 0.117] and [—1.3, 2.6] respectively. 

Keywords: cosmology:theory - chaplygin gas - dark energy - dark matter - CMBR 
- SNela 



1 Introduction 

It is now well established that the expansion of the universe is accelerating. The direct 
evidence for acceleration comes from the Hubble diagram of Type la Supernovae (SN la) pQ. 

*getsethi@physics.du.ac.in 



1 



The observational results of SN la together with the anisotropy of cosmic microwave back- 
ground radiation (CMBR) power spectrum and clustering estimates show that our universe 
is mainly made up of two components: dark matter and dark energy. The nature of dark 
matter and dark energy is not well understood. The dark matter contribute one-third of the 
total energy density of the universe. The dark energy, which is self interacting, unclustered 
fluid with large negative pressure, contributes roughly two-third of the total energy density 
of the universe (for latest review see [2]). 

The simplest and the most favored candidate for the dark energy is the cosmological 
constant (A). Consequently, several models with a relic cosmological constant (AC DM) 
have been used to describe the observed universe. However most of them suffer from severe 
fine tuning problem [21 d|- There are also models of dark energy with varying equation of state 
[I]. Another possibility for dark energy is quintessence which involve a slowly evolving and 
spatially homogeneous scalar field 0, or two coupled fields |7j. However , quintessence 
models also suffer from fine-tuning problem. This problem is usually highlighted as the "why 
now", that is, why does the dark energy start dominating over the matter content of the 
Universe recently. 

As an alternative to both the cosmological constant and quintessence, It is also possible 
to explain the acceleration of the universe by introducing a cosmic fluid component with an 
exotic equation of state, called Chaplygin gas [HUH]- The attractive feature of such models 
are that they can explain both dark energy and dark matter with a single component. The 
equation of state for the Chaplygin gas is P = —A/p, where A is a positive constant. A 
more generalized model of Chaplygin gas is characterized by an equation of state 

Pch = ~ (!) 

Pch 

where a is a constant in the range < a < 1 (the Chaplygin gas corresponds to a = 1). By 
inserting eq.(l) into the energy conservation law we get the expression for the energy density 
as PH 

^ = (^ + -^))^ (2) 

where a is the scale factor and B is a constant of integration. As can be seen from the 
above equation, the Chaplygin gas behaves like a non relativistic matter at early times while 
at late times the equation of state is dominated by a cosmological constant 87rGA 1 ^ 1+a ^ 
which thus leads to the observed accelerated expansion. The Chaplygin gas models have 
been found to be consistent with SNela data [TT] . CMB peak locations J2] and other ob- 
servational tests like gravitational lensing, cosmic age of old high redshift objects etc.[13j , 
as also with combination of some of them [T3]. It has been shown that this model can be 
accommodated within the standard structure formation scenarios (HlEBEi]- Therefore the 
Chaplygin gas model seems to be a good alternative to explain the accelerated expansion of 
the universe. However the Chaplygin gas model produces oscillations or exponential blowup 
of matter power spectrum that are inconsistent with observations [16, . 

Recently a variable Chaplygin gas model was proposed [T7] and constrained using SNela 
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"gold" data ^K|- I n this letter, we try to constrain the variable Chaplygin gas model pa- 
rameters using location of CMB peaks and SNela observations. We give the basic formalism 
of the model in section 2. Section 3 gives general formulae for CMB peaks. The best fit 
parameters have been calculated in section 4 by reduced \ 2 minimization of the distance 
modulus, followed by discussion and conclusions in section 5. 



2 Variable Chaplygin gas model 

We consider the proposed variable Chaplygin gas [Hj characterized by the equation of state: 

P+=-4& (3) 

Pch 

where A(a) = A a~ n , is a positive function of the cosmological scale factor a. A Q and n 
are constants. Using the energy conservation equation in flat FRW universe and eqlHJ the 
variable Chaplygin gas density evolves as: 



6 Aq B 

Pch(a) = \- ^ + - 4 

\ b — n a n a 

where B is a constant of integration. 

The original Chaplygin gas scenario is restored for n = and the gas behaves initially as dust- 
like matter and later as a cosmological constant. However, in the present case the Chaplygin 
gas evolves from dust dominated epoch to quintessence in present times (for details see ^7] 
and references therein). 

The Friedmann equation for a spatially flat universe reads as 

H 2 = — P (5) 

where H = a/a is the Hubble parameter. Therefore, the acceleration condition a > is 
equivalent to 

6-nJ V 

which requires n < 4. This gives the present value of energy density of the variable 
Chaplygin gas 



Pcho = \l — ^— A + B (7) 
6 — n 



where a = 1. Defining 



the energy density becomes 



1/2 



Pch(a) = Pch0 + (9) 



3 Location of CMBR peaks 

The location of peaks are very sensitive to the variations in the parameters of the model and 
hence serve as a sensitive probe to constrain the cosmological parameters and discriminate 
among various models [THl ED] - The peak locations are set by the acoustic scale I a, which 
can be interpreted as angle subtended by the sound horizon at the last scattering surface. 
This angle (say 9a) is given by ratio of sound horizon to angular diameter distance, cIa of 
the last scattering surface: 

e A = a M^M (10) 

where c s is the speed of sound given in plasma by c s = 1/(3(1 + R)) 1 ' 2 and R = 3pb/ p 1 
corresponds to the ratio of baryon to photon density. For the flat (k = 0) FRW model, the 
Acoustic scale I a = t^/Oa is 

, To _ T ls n , \ 

l A = 7T (11) 

C s T is 

where r and r is are conformal time today and at the last scattering surface. 

For our calculations we assume: present value of scale factor a = 1, the value of scale 
factor at the last scattering surface ai s = 1/1100, h = 0.71, the density parameter for 
radiation and baryons at present f2 7o = 9.89 x 10~ 5 , Qbo = 0.05, average sound speed 
c s = 0.52 and spectral index for initial energy density perturbations, n — 1. 

The location of the i th peak in the spectrum is given by 

U = U{i ~ Si) (12) 

where the phase shift <5j, caused by the plasma driving effect, is solely determined by the 
pre-recombination physics. We can approximate it with its value from standard cosmology 

Sl ~ .267[^^] ai (13) 

where r(zi s ) = p r (zi s ) / p m (zi s ) = fi r o(l + Z[ s )/Q m0 is the ratio of radiation and matter 
density at the decoupling epoch. For the third peak, 



5-, w 0.35 



r{zi s ) 
L 0.3 



o.i 



(14) 



4 



The Friedmann equation, using Eqn. El for a variable Chaplygin gas becomes 



H 2 



8nG 



p r0 (i + zf + Pb0 {i + zf + Pch0 n m (i + zf + (i - a m )(i + z ) n 



1/2 



(15) 



where p r0 and pbo are the present values of energy densities of radiation and baryons respec- 
tively. 



Using 



and 



Eqn. ITB1 becomes 



PrO 


f2 r o 






PchO 


QchO 


1 


— fi r o — fifeo 


PbO 


^bO 




QbO 


PchO 


&<chO 


1 


— fi r0 — fifeo 



where 



X 2 (a) 



n 



rO 



H = ^ cm H 2 a^X 2 (a) 

QbO a 



, O 1 _ O \ 1/2 

+ a —r- H 
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Using the fact that H 2 = C ~ A '^ 2 



we get 



dr 



da 



nl^H X(a) 



so that Eqnll II becomes 



71 



1 da 



X(a)\J X{a) 



da \-i 



(16) 
(17) 



(19) 



(20) 



(21) 



The ratio of radiation to matter density at the last scattering surface, ri s in Eqnll3l becomes 



af s ( 1 - Qro - ^ 



+ 



(l-O m ) \ 
a n J 



1/2" 



(22) 



4 x 2 minimization 

We have used the SNe la data to constrain the parameters of the variable Chaplygin gas 
model. Here, to find out the luminosity distance we have taken into account the contributions 
from radiation and baryons together with the Chaplygin gas. 
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Using the Friedmann equation ^3], in flat universe, the luminosity distance can be ex- 
pressed as 



di 



da 



and the distance modulus 



^ = 51og^ + 42.38 
cn 

The best fit parameters are determined by minimizing 



X 



fhh Pobs 



where 



at 



i 

C -E4 



- 7 ^-fCi + -lnl0) -2\nh 

G2 V 5 / 

Wth ~ Vobsl 



(23) 



(24) 



(25) 



(26) 



(27) 



We obtain the minimum \ 2 = 174.36 for Q m = 0.22 and n = —2.8 which are quite close 
to those obtained by [T2j. The contour corresponding to 99% CL in figure 1 & 2 gives the 
range of Q m = [0.0,0.36] and n = [—41.3,2.8] as tabulated in Table 1. 



5 Discussion 



The Variable Chaplygin gas behaves like non-relativistic matter during the early times and 
has the ability to drive the universe into an accelerated expansion in the recent times. From 
SNe la Gold data, this model is consistent with wide range of values for parameters Q m and 
n. The best fit of the model obtained bv Quo fc Zhang [T7j gives Q m = 0.25 and n = —2.9 
with xLin = 173.88. They performed the Xmin ^ es * whh the gold data set [T8j using 



H 2 



8ttG 



PchO 



fi m (l + ^) 6 +(l-fi m )(l + ^) 



1/2 



(28) 



We repeated the test using the same SNe la Gold data set but took into account contribution 
from matter as well as radiation, as given in Eq. The results obtained by us give 

Xmin = 174.36 at Q m = 0.22 and n = — 2.8. We find that matter and radiation do 
not contribute significantly to the expansion rate at redshifts for which the supernova data 
is available but we include them for completeness of our analysis. At 99% confidence level 
the Xmin i s obtained for values of parameters Q m and n lying in the range [0.0, 0.36] and 
[-41.3, 2.8] respectively. 
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The Figures 1 and 2 show the la, 2a and 3a x 2 contours of the SNe la Gold data and 
the location of CMBR peaks (using results of WMAP) at la level in the Variable Chaplygin 
Gas model. In Figure 1, the contours 1 and 2 represent, for the first peak, the la lower 
(h = 291.3) and the upper (h = 220.9) bounds on I obtained from WMAP data. In 
conclusion, the la contours of the first peak of the WMAP combined with SNe la Gold 
sample restrains the parameter space to Q m = [0.017, 0.117] and n = [—1.3, 2.6]. In 
Figure 2, the contours 1 and 2 represent, for the third peak, the la lower (7 3 = 802.0) and 
the upper (Z 3 = 816.0) bounds on I obtained from WMAP data. The la contours of the 
third peak of the WMAP combined with SNe la Gold sample restricts the parameter space 
to Q m = [0.018, 0.091] & n = [—0.2, 2.8]. The results have been tabulated in Table 1. 

In Figures 3 and 4 , we give the variation of the location of first and third peak respec- 
tively with n for different values of Q m . The horizontal lines in these two figures show the 
observational la bounds for the two peaks. The bound on Q m from the first peak is Q m < 0.2 
and from the third peak is fl m < 0.18. 

Recently Guo and Zhang [22] have put constraints on the Variable Chaplygin Gas using 
X-Ray gas mass fractions in galaxy clusters and SNe la Gold data set. They have shown that 
the consistent range of Q m and n are [0.043, 0.068] and [1.18, 2.03] (la error bar) respectively. 
Our results based on location of WMAP 3 rd peak and Snela at the 3a level indicate the 
range [0.018, 0.091] and [—0.2, 2.8] for Q m and n respectively which accommodates the 
results obtained from X-Ray gas mass fractions in galaxy clusters and SNe la [22] • Though 
the original Chaplygin Gas is observationally ruled out, the Variable Chaplygin Gas can be 
further explored. 



Table 1. Bounds on Q m & n. 



Method 


Reference 




n 


SNeIa(3<r) 


This paper 


22 +0 14 


O q+5.6 
^•°-38.5 


WMAP 1 st Peak 
+ SNeIa(3<r) 


This paper 


[0.017, 0.117] 


[-1.3, 2.6] 


WMAP 3 rd Peak 
+ SNeIa(3a) 


This paper 


[0.018, 0.091] 


[-0.2, 2.8] 
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Figure 1: The combined WMAP li and SNela parameter space 
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Figure 2: The combined WMAP l 3 and SNela parameter space 
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Figure 3: Variation of li with n for different Q m . The horizontal lines correspond to the 
observational la bound on l\ 
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Figure 4: Variation of l 3 with n for different Q m . The horizontal lines correspond to the 
observational la bound on l 3 
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